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Abstract 

Background: Syn6scXy\y type 1 (SDI) is an autosomal dominant limb malformation characterized in its classical form by 
complete or partial webbing between the third and fourth fingers and/or the second and third toes. Its four subtypes (a, b, 
c, and d) are defined based on variable phenotypes, but the responsible gene is yet to be identified. SDI -a has been 
mapped to chromosome 3p21.31 and SDl-b to 2q34-q36. SDl-c and SDl-d are very rare and, to our knowledge, no gene 
loci have been identified. 

Methods and Results: \n two Chinese families with SDl-c, linkage and haplotype analyses mapped the disease locus to 
2q31-2q32. Copy number variation (CNV) analysis, using array-based comparative genomic hybridization (array CGH), 
excluded the possibility of microdeletion or microduplication. Sequence analyses of related syndactyly genes in this region 
identified c.917G>A (p.R306Q) in the homeodomain of H0XD13 in family A. Analysis on family B identified the mutation 
C.9160G (p.R306G) and therefore confirmed the genetic homogeneity. Luciferase assays indicated that these two 
mutations affected the transcriptional activation ability of H0XD13. The spectrum of H0XD13 mutations suggested a close 
genotype-phenotype correlation between the different types of HOXD /3-Syndactyly. Overlaps of the various phenotypes 
were found both among and within families carrying the H0XD13 mutation. 

Conclusions: Mutations (p.R306Q and p.R306G) in the homeodomain of H0XD13 cause SDl-c. There are affinities between 
SDl-c and synpolydactyly. Different limb malformations due to distinct classes of H0XD13 mutations should be considered 
as a continuum of phenotypes and further classification of syndactyly should be done based on phenotype and genotype. 
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Introduction 

Syndactyly (SD) is a digital malformation in which adjacent 
fingers and/or toes are webbed because they fail to separate 
during limb development. It is one of the most common hereditary 
limb malformations with a prevalence of 3-10 in 10 000 births, 
though higher estimates ranging from 10-40/10 000 have been 
reported [1-2]. The current classification scheme defines nine 
types syndactyly on the basis of phenotype and genotype, which is 
helpful in the understanding of syndactyly malformation and in 
defining its affmity witii other digit anomalies, including Polydac- 
tyly, oligodactyly, and brachydactyly [3]. Of all the known non- 
syndromic syndactylies, SDI is one of the most common types and 
follows an autosomal dominant inheritance pattern with variable 
phenotypes. Webbing may affect fingers and/ or toes, be unilateral 
or bilateral, cutaneous or bony, and reach the level of the nail or 
solely affect the proximal segments of the digits. Based on clinical 
observations and current classification, SD 1 can be further divided 
into four subtypes (a, b, c, and d). 



SDI -a (l\'eidenreich type; zygodactyly; OMIM 609815) was 
mapped to chromosome 3p21.31 and characterized by bilateral 
webbing of 2/3 toes without hand anomalies. The features of SDl- 
b (Lueken type; OMIM 185900) are bilateral cutaneous 3/4 
fingers and webbing of 2/3 toes. It was mapped to chromosome 
2q34-q36 and the chromosome locus was designated as the SDI 
locus [4-5]. SDl-c, also named Montagu type or 3/4 fingers 
syndactyly, is characterized by unilateral or bilateral cutaneous or 
bony webbing of 3/4 fingers and occasionally of 4/5 fingers. It is 
recognized that the feet are not involved in this subtype [6-7]. 
SD 1 -d (CastiUa type) has only been observed in an epidemiological 
study [2] and is characterized by bilateral cutaneous webbing of 4/ 
5 toes. To date, littie is known about these four subtypes and the 
gene locus for the rare SDl-c and SDl-d subtypes has not been 
identified. 

HOX genes encode a set of highly conserved transcription 
factors that control cell fate and the regional identities along the 
primary body and limb axes in metazoans [8]. There are 39 HOX 
genes arranged in four separate clusters (HOXA-D) in mammalian 
genomes. HOX proteins bind specific DNA sequences via their 
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homeodomains, and are thought to regulate overlapping sets of 
target genes, although the molecular pathways controlled by HOX 
proteins are only starting to be unveiled [9] . II0XD13, the gene at 
the 5' end of the HOXD cluster, is the first //OX gene known to be 
linked to human developmental disorders. It has two coding exons: 
exon 1 with the impc-rft-ct GCN (where N denotes A, C, G, or T) 
triplet repeats encoding the N-terminal region with a 15-residue 
polyalanine tract, and exon 2, which contains the homeobox 
region encoding the C-terminal portion with a 60-residue 
homeodomain [10]. Mutations in the human H0XD13 gene 
could lead to diversiform phenotypes, including synpolydactyly 
type 1 (SPDl; OMIM 186000), brachydactyly types D (BDD; 
OMIM 113200) and E (BDE; OMIM 113300), SD5 (OMIM 
186300), a novel brachydactyly-syndactyly syndrome (BD-SD; 
OMIM 610713), or VACTERL association (OMIM 192350). 

In this study, our aim was to further characterize the SDl-c 
phenotype based on two Chinese families with 3/4 fingers 
syndactyly, ^\'e discovered that two missense mutations in codon 
306 in the homeodomain oi HOXD 13 underKe SDl-c. 

Materials and Methods 

Ethics Statement 

The study was approved by the Third Military Medical 
University institutional review board. A written informed consent 
was obtained from all participants and/or guardians on the behalf 
of the minors/ children participants. 

Families 

Two Chinese Han families (Figure lA and IB) were ascertained 
in the present study. After physical examination, digital photo- 
graphs and radiographs were taken to document the limb 
phenotypes in some of the individuals. 

Family A (FA) consisted of 74 individuals including 33 affected 
individuals (12 females and 21 males) of five generations. The 
proband (FA-IV26) had complete bilateral webbing of 3/4 fingers, 
where the webbing included a fusion of nails according the 
phenotype and radiographs. His feet were normal and no other 
abnormalities were noted. The phenotypes varied from partial 
unilateral to complete bilateral c:utaneous webbing of 3/4 fingers 
in the other affected individuals. Among the 18 affected 
individuals available for phenotypic evaluation, 13 had the same 
phenotype as the proband, 2 had unilateral cutaneous webbing of 
3/4 fingers (FA-IV4, FA-V2) and only FA-V9 had a mild SPD 
phcnot)'pe of the hands with normal feet. Notably, two affected 
individuals (FA-II13, FA-1II22) presented with an atypical abnor- 
malit}' of the feet but without hand anomalies. FA-1II3 had a 
unilateral adduction deformity of the toes. FA-III22 had unilateral 
partial cutaneous webbing of 2/3 proximal toes (Figure lA, 
Figure 2, Figure SI and Table 1). 

Family B (FB) consisted of 13 members, including 7 affected 
individuals females and 2 males). All members were available for 
clinical evaluations. Five affected individuals had a complete 
bilateral webbing of 3/ 4 fingers with normal fe(;t. One of them 
(FB-IIIl) also had clinodactyly of the fifth finger. In addition, FB- 
III4 presented unilateral cutaneous webbing of 2/3 fingers with 
normal feet. Notably, FB-IIL5 had a unilateral synpolydactyly of 
4/5 toes and unilateral cutaneous webbing of 3/4 fingers 
(Figure IB, Figure 2, Figure SI and Table 1). 

Linkage and Haplotype Analyses 

As we known, linkage analysis remains a very good choice for 
large family and could clearly reflect the disease locus. Short 
tandem repeat (STR) typing was performed firstly in present 



families. Blood samples were collected and genomic DNA was 
extracted from 29 and 1 1 members of family A and family B, 
respectively. We carried out a search of the whole chromosome 2 
(30 microsatellite markers) and chromosome 3p21.31-3p22.31 
region (2 microsatellite markers) on the base of Linkage Mapping 
Set v2.5. Three additional microsatellites (D2S2981, D2S2314, 
and D2S324) were selected using the UCSC Genome Browser on 
Human Mar. 2009 Assembly (http://www.genome.ucsc.edu) to 
further refine the critical region of the disease locus. Two-point 
LOD scores were calculated using the MLINK program of the 
FASTLINK package, assuming that the disease in the family is 
inherited in an autosomal dominant mode with penetrance of 
0.95, the disease-aUele frequency is 0.0003, and allele frequencies 
are equal at aU the marker loci. 

Microarray-based Comparative Genomic Hybridization 

Array-based comparative genomic hybridization (array CGH) 
was carried out using the Roche NimbleGen Genome-Wide array 
CGH 3*720K containing over 720,000 copy-number probes. 
Genomic DNA samples were genotyped at the CapitalBio 
Corporation (Beijing, China) with the CGH array in accordance 
with the manufacturer's protocols. Genotype calling, genotyping 
quality control, and copy number variation (CNV) identification 
were performed with the Roche NimbleGen SignalMap software. 
The data accessed at Gene Expression Omnibus (GEO, No: 
GSE55181). 

Mutation Screening 

Exons and splice sites of syndactyly-related genes H0XD13, 
GJAl and zone of polarizing activity regulatory sequence {ZRS) 
were amplified by polymerase chain reaction (PGR) using gene- 
specific primers. Amphfied fragments were sequenced using an 
ABI 3130 automated sequencer (Applied Biosystems, Foster City, 
CA, USA). 

Review of the H0XD13 iVlutations 

The name of H0XD13 mutations in related studies is confusing 
because it is uncertain whether Met-1 or Met-9 is the initiator at 
present. The primitive sequence AAC5 1635.1 contains 335 amino 
acids and the newly accession NP_000514.2 contains 343 amino 
acids, which could result from erroneous and uncertain initiation. 
Thus, we reviewed the H0XD13 mutations based on sequence 
NP_000514.2 to better understand its mutation spectrum. 

Functional Significance Prediction 

Three online programs: PolyPhen2 (http://genetics.bwh. 
harvard.edu/pph2), SIFT (http://sift.jcvi.org), and Swiss-Model 
(http:// swissmodel.expasy.org), were used to predict the functional 
significance of mutations. 

Piasmld Construction and Luciferase Assays 

The human EPHA7 gene promoter of 660 bp (from —580 to 
+80), which contains a HOXD 13 binding site, was obtained by 
PGR from human genomic DNA and insert into the Xhol and 
Hindlll sites of a pGL3-basic vector (Promega, Madison, WI, 
USA) to generate a pGL3-EPHA7 reporter construct. Human 
HOXD 13 of pReceiver-M02 expression clone was obtained from 
GeneCopoeia (RockviUe, MD, USA). Site-directed mutagenesis 
was performed with appropriate primers to generate HOXD 13 
carrying die C.9160G (p.R306G), c.917G>A (p.R306Q) or 
c.964A>C (p.I322L) mutations using the QuickChange Lightning 
Site-DirectedMutagenesis kit (Strata-gene, La Jolla, CA, USA). 
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Figure 1. Pedigree, Kiaplotype analysis, and mutation analysis of the families in this study. A: Pedigree and disease-haplotype 
segregation of family A. Blackened symbols represent affected Individuals with an abnormal limb phenotype. White symbols represent Individuals 
with a normal limb phenotype. Circles and squares Indicate females and males, respectively. The arrows Identify the proband and the disease- 
haplotype Is boxed. B: Pedigree and disease-haplotype segregation of family B. C: H0XD13 missense mutations In the two families. 
dol:1 0.1 371/journal.pone.00961 92.g001 



Presence of the desired base changes was verified by DNA 
sequencing. 

The 293FT cells were cultured in 1640 supplemented with 10% 
FBS, and seeded in 24-weU tissue culture plates 24 h prior to 
transfection at about 80% confluence. Plasmids were transfected 
with lipofectamine2000 (Invitrogen) according to the manufactur- 
er's instructions. HOXD13 expression vectors or the control 
vectors were cotransfected with pGL3 reporter vectors along with 
pRL-TK vectors (ReniUa luciferase, Promega) using lipofecta- 
mine2000. Cell was collected 24 h after transfection and luciferase 



activities were measured using Dual-Luciferase Reporter Assay 
System (promega). Activity was defined as Firefly/ReniUa ratio. 

Results 

Segregation Analysis 

In family A, the syndactyly did not skip any generations, 
indicating that the disease is of autosomal dominant inheritance. 
Segregation analysis of the pedigree confirmed this mode of 
inheritance. Results from the two-point linkage analysis provided 
strong evidence that the SD 1 -c disease locus in this family linked to 
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STR markers of chromosome 2q31-2q32. The maximum LOD 
score calculated was 4.88 at a recombination fraction of 0.00 and 
penetrance of 0.95 at the D2S335 locus. The LOD scores of 
2q33.3-q36.3 and 3p21.31-3p22.1 regions were insufficient to 
support the hnkage relationship. Haplotype analysis demonstrated 
that all of the 17 affected members examined had a common 
haplotype of 3-4-1-2 for four STR marker loci: D2S335, 
D2S2981, D2S2314, and D2S324. Based on haplotype analysis, 
recombination events were observed in FA-IV 1 1 and FA-V 1 . One 
unaffected individual (FA-Vl) carried a common disease-associat- 
ed haplotype and the sequencing results demonstrated that the 
individual carried a mutation. (Figure lA and Table 2). 

In family B, the common region was also confirmed. A 
haplotype of 5-5-6-1 of this four STR marker loci (D2S335, 
D2S2981, D2S2314 and D2S324) was found in aU patients. The 
maximum LOD score in family B was 1.79. (FigurelB and 
Table 2). 

CNV Analysis 

To determine whether the SDl-c phenotype was caused by an 
unknown microdeletion or microduplication, we performed a 
genome-wide high-resolution CNV scan of an affected individual 
(male, FA-IV4) and a healthy control (male FA-IV8). However, we 
did not detect any potential pathogenic CNVs in the critical 
regions (GSE55181), including the ^i?.? region which is the cause 
of SD4. 

Mutation Analysis 

Direct sequencing of H0XD13 revealed a heterozygous G-to-A 
transition in exon 2 at position 9 1 7 of the coding sequence in the 
proband of family A. This base change converted amino acid 306 
from arginine to glutamine. The same base change was identified 
by direct sequencing in other affected family members and one 
unaffected individual (FA-Vl). In the proband of family B, there 
was a G-to-G transversion at position 916, thus leading to the 
change of amino acid 306 from arginine to glycine. These two 
missense mutations were not found in unaffected family members 
except FA-Vl or in the 100 unaffected healthy control subjects 
(Figure IG). No mutations were found in the GJAl gene and ^RS 
region. 

A Spectrum of H0XD13 Mutations 

Our results provide new information and enlarge the spectrum 
of known H0XD13 mutations. To date, more than 20 different 
variants of the H0XD13 gene have been reported from different 
families [9,11-28]. Some mutations with unspecified phenotypes 
have also been found in sporadic cases [29] . The present cases in 
combination with those studied previously suggest that mutations 
in the polyalanine tract and homeodomain are in the majority of 
families with H0XD13 mutations. These data suggest that 
mutations in these regions of the protein underlie the majority of 
//OXZ)75-associated cases and that this region may be particularly 
important for protein function (Figure 3A and Table 3). 

Functional Analysis 

The p.R306Q_and p.R306G mutations altered the amino acids 
at residue 3 1 of the homeodomain. An alignment of HOXD 1 3 
protein sequences showed that this position is highly conserved 
among many different species (Figure 3B). Thus, this amino acid 
appears to play an important role in the structure and function of 
the HOXD 13 protein. 

The R306Q_ and R306G substitutions were rated as 'probably 
damaging' by PolyPhen2 analysis, with a score of 1 .00 (sensitivity: 
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Figure 2. Photographs and radiographs of individuals. Phenotypic spectrum of the two families showing typical SDI-c bilateral 3/4 fingers 
syndactyly (A, B, H, I and J) and unilateral 3/4 fingers syndactyly (G, iVl and N); transformative SDI-c bilateral 3/4 finger syndactyly with right fifth 
finger clinodactyly (K and L); mild SPD, bilateral 3/4 fingers syndactyly with phalanx duplication (E and F) and unilateral 3/4 fingers syndactyly with 
4/5 fingers synpolydactyly (O and P). Other phenotypes included isolated unilateral partial webbing of 2/3 toes like SD1 -a (D) and isolated unilateral 
adduction deformity of toes (C). 
doi:1 0.1 371 /journal.pone.00961 92.g002 



0.00; specificity: 1 .00). Other misssense mutations were also rated 
as 'probably damaging'. However, Swiss-Model analysis indicated 
that the overall structures of the mutant proteins did not differ 
substantially from the wild-type protein. 

Luciferase assays were conducted to determine whether the 
mutation affected the capability of HOXD 1 3 protein to activate 
transcription. The pGL3-EPHA7 luciferase reporter constructs 
were tested and the c.964A>C (p.I322L) mutant that converts 
residue 47 of the homeodomain from isoleucine to leucine was also 
examined (approx. 58% of wild type). The C.9160G (p.R306G) 
and C.9170A (p.R306Q) mutants also showed significantly 
diminished stimulation compared with the wild-type control 
(approx. 62% and 60% of wild-type, respectively). Thus, these 
two mutations affected the capacity of HOXD 13 to activate 
transcription (Figure 4). 

Discussion 

SDl-c was first described by Motagia in 1953 and is 
characterized by unilateral or bilateral, cutaneous or bony fusion 
of 3/4 fingers with normal feet [7]. Another large Chinese family 
with SDl-c, reported by Hsu, had 23 affected subjects demon- 



strating variable degrees of bilateral osseous fusion of 3/4 or 3/ 4/ 5 
fingers. Only one subject had partial fusion of the toes from the 
third to the fifth [6] . The disease followed an autosomal dominant 
inheritance pattern with phenotypic heterogeneity and incomplete 
penetrance in the families of the present study. Most affected 
individuals exhibited isolated bilateral syndactyly of 3/4 fingers, 
while a small number were unilateral (FA-IV4, FA-V2, FB-III4). 
Other particular phenotypes were also observed as follows: (1) 
unilateral abnormal foot without hand anomalies; FA-III3 only 
presented with unilateral adduction malformation of the toes and 
FA-III22 had a unilateral partial webbing of 2/3 toes overlapping 
the typical appearance of SDl-a. (2) Occasional mild SPD 
phenotype; FA-V9 had bilateral SPD (3/4 fingers) without foot 
anomalies; FB-III5 had unilateral 3/4 fingers syndactyly with a 
mild unilateral phenotype of SPD (4/5 toes). (3) A rare individual 
with forme fruste (FA-Vl) was observed. The families in our study 
confirmed the genetic homogeneity of previous SD 1 -c studies and 
the typical SDl-c phenotype was observed in almost all affected 
individuals with only a small number of affected individuals 
exhibiting an alternative appearance. 
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Figure 3. Schematic presentation of the HOXD1 3 structure and sequences alignment analysis. A: Schematic presentation of the I-I0XD13 
structure and annotated mutations identified in families. B: Partial amino acid sequences alignment in homeodomain of H0XD13 among several 
species. The position of mutated amino acid of our cases is indicated by the black box. 
doi:1 0.1 371/journal.pone.00961 92.g003 



There were some common associations between SDl and SPD. 
The typical feature of SDl-c is isolated syndactyly of 3/4 fmgers, 
while the cardinal features of SPD are the webbing of 3/4 fingers 
and 4/5 toes, with partial or complete digital duplication within 
the syndactylous web. The webbing phenotype (3/4 fmgers) was 
recognized as a milder SPD phenotype in several SPD families. In 
addition, the SDl -a phenotype (isolated webbing of 2/3 toes) was 
occasionally found in SPD families [32] [33] . The close association 
between SDl-c and SPD inspired us to determine the disease- 
associated gene for SDl-c. The evidence provided by linkage and 
haplotype analyses revealed that affected members of family A had 
a common haplotype of 3-4-1-2 for four STR marker loci, while 
family B had 5-5-6-1. The different haplotypes suggested that 
these two families were derived from different ancestors; therefore, 
the region was localized to 2q31-2q32. The evidence from 
overlapping phenotypic and genetic analyses indicated that 
abnormalities in H0XD13 may be the common cause of SDl-c 




Bnpty WT R31Q R31G I47L 



pGLS-BSHA/ 

Figure 4. Transcriptional activity of wild-type and mutant 
HOXD13 proteins at the human EPH7A promoter. The R31Q 
(p.R306Q), R31G (p.R306G) and I47L (p.l322L) mutants show a 
significantly diminished stimulation compared with the wild-type 
control at approximately 62%, 60% and 58% of the wild-type H0XD13, 
respectively. Bars represent firefly/Renilla luciferase ratios for the 
different constructs. 
doi:1 0.1 371/journal.pone.00961 92.g004 



and SPD. Finally, the result of the direct sequencing of H0XD13 
confirmed our hypothesis. Two H0XB13 mutations, p.R306Q^ 
and p.R306G, were identified in the families of the present study, 
respectively. At the same time, CNV analysis using array 
CGH excluded the cause of an vmknown microdeletion or 
microduplication, such as in the region. In our previous 

study of SD4, one individual that carried a ^RS duplication 
exhibited a postaxial Polydactyly of the right foot similar to SPD 
[34]. 

The H0XD13 gene is associated with SD-lc in the presented 
families. Other limb malformations associated with mutations in 
H0XD13 included SPDl, BDD, BDE, SD5, BD-SD syndrome, 
and VACTERL. They provide a large dataset for analysis of the 
most common mutations and of the different types of mutations 
found, thus generating a preliminary genotype-phenotype corre- 
lation. Our present report and the observations of others are an 
indication of the existence of genotype-phenotype correlation for 
the limb morphopathies caused by H0XD13 mutations, especially 
in the homeodomain. 

The homeodomains of the HOX genes are DNA-binding motifs 
composed of a flexible N-terminal region followed by three a- 
helices and are highly conserved among different species of the 
HOXD13 protein. In previous studies, five missense mutations 
have been identified in the homeodomain (R298W, R298Qj 
S308C, I314L and Q317R based on AAC51635.1, also known as 
residues R31W, R31Q, S41C, I47L and Q50R of the homeodo- 
main), each producing distinctive phenotypes. Interestingly, the 
two mutations (R306Q and R306G based on NP_000514.2, also 
named R31W and R31G) located in the same codon (306) of 
present two families also exhibit a different phenotype from that 
reported previously. 

The R31W mutation was previously in one Danish family 
reported by Debeer et al. Fifteen affected individuals had bilateral 
fifth finger clinodactyly. Three of those also had mild unilateral 
SPD, which was radiologicaUy documented in two cases as 
cutaneous syndactyly of 3/4 fmgers, with a duplicated distal 
phalanx in the syndactylous web. In addition, one subject was 
clinically unaffected but an obligate mutation carrier [21]. The 
R3 1 Q_ mutation was documented in a Chinese family by Wang et 
al. In that study, three individuals presented with an isolated 
bilateral cutaneous webbing of 3/4 fingers. Three had bilateral 
clinodactyly of the index fingers and one had bilateral clinodactyly 
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of the second toe. In the family A of our present study, the RSIQ 
mutation was observed, where twenty-three affected individuals 
exhibited isolated bilateral or unilateral 3/4 fingers syndactyly. 
Three affected individuals exhibited appearance of mild SPDl, 
SDl-a or isolated unilateral toe adduction. The R31G mutation 
was observed in family B and six affected individuals had isolated 
bilatersd or unilateral 3/4 syndactyly, one affected individual had 
unilateral 3/4 fingers syndactyly with unilateral mUd phenotype of 
SPD. Moreover, carriers of the S41C and I47L mutations 
displayed a unique brachydact>4y-polydactyly syndrome, overlap- 
ping with BDE, with considerable variations in severity and 
expressivity. Carriers of Q50R displayed another special pheno- 
type of SD5. 

This phenotypic distinction might arise from the specific 
impact of different amino acids on protein stability or its 
affinity for binding DNA recognition elements. An R31P 
missense mutation in the homeodomain of human MSXl led 
to partial or complete loss of function because of a perturba- 
tion of the structure and diminished ability compared with the 
wild-type gene [35]. EPHA7 gene is downstream of HOXD13 
and acts during limb development. HOXD13 can bind to a 
single evolutionarily conserved site within the EPHA7 promot- 
er and activates its transcription [36]. Transactivation assays 
also demonstrated that some mutations affected the transcrip- 
tional activation ability of HOXD13. Both the R31Q,and I47L 
mutants showed a significantly reduced transcription activation 
ability of the EPHA7 promoter compared with the wild type 
[11]. The Q50R mutation resulted in a more reduced 
activation than the I47L mutation [24]. The proximal 
cartilages in embryos with ectopic expression of HOXD 1 3 
(IQ_N), an artificial mutant carrying alanine substitutions at 
positions 47(1), 50(Q) and 51(N) of the homeodomain, which 
renders it unable to bind DNA, were shorter than in the 
embryos with expression of wild-type HOXD 13 [37]. Our 
Luciferase assays demonstrated that R31 and G31 mutants also 
changed the DNA-binding ability. All missense mutations in 
HOXD 13 were rated as 'probably damaging' by PolyPhen2 
analysis. This suggests that mutations in the homeodomain 
result in a partial loss of function in the transcription- 
regulating functions of HOXD 13. Therefore, an alteration of 
arginine would affect the ability of HOXD 1 3 to make the salt 
bridges normal!)' formed by R3 1 . Tliis plienotypic distinction 
might arise from the impact of different amino acid residues on 
the stability of the protein structure or the DNA-binding ability 
of the homeodomain to the DNA recognition elements. 
Alternatively, some unidentified genetic factors, except the 
impact of CNV, might correlate with H0XD13 mutations 
resulting in varied phenotypes. Therefore, the substitution 
would affect the structure and lead to loss of function. 

We re-examined the proposed genotype-phenotype correla- 
tion for HOXD 13 (Figure 3 and Table 3). The polyalanine 
repeat expansions and frameshift deletion mutations usually 
give rise to typical SPD phenotypes, and the homeodomain 
missense mutation cluster for brachydactyly types [38]. Close 
genotype-phenotype correlations were also observed between 
different types of //OXZ)75-Syndactyly. Overlap of various 
phenotypes could be found among and within families carrying 
H0XD13 mutation. The same minor malformations may be 
associated with different clinical abnormalities of the limb. For 
example, soft-tissue syndactyly of 2/3 toes were observed in 
famihes with mild SPD, BD-SD syndrome, SD5 [24], SDl-a, 
SDl-b or SDl-c. Additionally, cutaneous syndactyly of 3/4 
fingers without digit duplication was observed in mild SPD, 



SD5, SDl-b or SDl-c. In the same codon of the same gene, the 
phenotype of R31 mutations in HOXD 13 exhibited an 
interesting continuous phenotype: from unilateral to bilateral, 
from SDl-c phenotype to isolated 3/4 finger SPD and typical 
SPD (3/4 finger SPD and 4/5 toe SPD), from mild SPD to 
severe SPD and other phenotypes. No definite line could be 
drawn between the milder SPD variants and the SDl 
phenotypes. Therefore, different limb malformations due to 
distinct classes of HOXD 13 mutations should be considered as 
a continuum of phenotypes. "HOXD 13 limb morphopathies" 
appears to be the most appropriate term to cover all limb 
malformations due to HOXD 13 mutations in accordance with 
Zhao et al., Malik et al., and Brions et al. [24-25,38]. 

Typing of SDl-c may be quite puzzling since phenotypic 
heterogeneity may lead to a clinical misclassification. Intra- 
familial and individual-specific clinical variability in SDl-c 
families may exhibit phenotypic overlap with other entities 
such as SDl-a, SPD or atypical limb malformations (e.g. 
clinodactyly). The same phenotype could be caused by 
different gene mutations and mutations in the same gene 
could also lead to different phenotypes such as that of the 
H0XD13 mutation. In addition, different missense mutations 
in the same codon exhibited phenotypic heterogeneity. 
Individuals with rare forme fruste further reflect the complex- 
ity of phenotype. The heterogeneity may partly reflect 
individual lifestyles and specific types of mutation. In addition, 
environmental factors and other genetic factors such as genetic 
polymorphisms linked to SD-associated genes and epigenetic 
regulation of the expressicm of these genes may be involved in 
these syndromes. Whole cxome sequencing of these families 
could provide important insights in the future. At present, 
more weight should be given to the phenotypes seen in large 
families for the correct typing of SD (and other developmental 
anomalies) and the most common phenotype in these large 
families should be considered firstly for typing [3] . 

In conclusion, mutations (p.R306Q_ and p.R306G) in the 
homeodomain of HOXD 13 could cause SDl-c. There are affinitis 
between the SD-lc and SPD phenotypes. Our study provides new 
information regarding the involvement of H0XD13 gene muta- 
tions in the clinical and genetic diversity of SD, and advances our 
understanding of human limb development. Different limb 
malformations due to the distinct classes of H0XD13 mutation 
should be considered as a continuum of phenotypes and further 
classification of syndactyly should be done based on phenotype 
and genotype. The term "HOXD 13 limb morphopathies" seems 
appropriate to cover all limb malformations due to H0XD13 
mutations. 
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